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<N : This is the summary of a week of very informative presentations on new ways 

to probe the Universe using gravitational detectors, space and ground based gamma 
■ ray telescopes, EeV air shower detectors and neutrino telescopes. Gamma ray bursts 

and active galaxies were hot theoretical themes in the multi-wavelength discussions. 
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1. Introduction 



Efforts are underway to qualitatively improve the instruments that can push astron- 
omy beyond GeV photon energy, to wavelengths smaller than 10 _14 cm, and map the 
sky in neutrinos and EeV cosmic rays as well as gamma rays. New gravitational 
wave detectors will explore wavelengths much larger than those of radio astronomy. 
The diffuse cosmic photon flux is shown in Fig. 1 as a function of photon energy 
and wavelength. The instruments which have collected the data shown in the figure 
span 60 octaves in photon frequency, from 10 4 cm radio-waves to 10~ 14 cm photons of 
GeV energy. This is an amazing expansion of the power of our eyes which scan the 
sky over less than a single octave just above 10 _5 cm wavelength. The new astron- 
omy, discussed at this conference, probes the Universe at new wavelengths, smaller 
than 10~ 14 cm and larger than 1 kilometer. Besides gamma rays, gravitational waves, 
neutrinos and the very high energy protons which are only weakly deflected by the 
magnetic field of our own galaxy, become the messengers from the Universe. As 
exemplified time and again, the development of novel ways of looking into space in- 
variably results in the discovery of unanticipated phenomena. As is the case with new 
accelerators, observing the predicted will be slightly disappointing. 

*Talk given at The XXXIInd Rencontre de Moriond, "Very High Energy Phenomena in the 
Universe", Les Arcs, France (1997). 
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Figure 1: Flux of gamma rays as a function of wavelength and photon energy. In 
the TeV-EeV energy range the anticipated fluxes are dwarfed by the cosmic ray flux 
which is also shown in the figure. 



The "New Astronomy" is sketched in Fig. 1 and covers: 

• gravitational waves with the commissioning of pairs of resonant bar detectors 
and LIGO, 

• gamma ray astronomy, with new satellite-borne detectors expanding from the 
GeV, into the hundreds of GeV energy region, while, at the same time, second- 
generation ground-based air Cherenkov telescopes reach down to thresholds of 
tens of GeV, 

• cosmic rays which may point at their sources with degree accuracy when their 
energy exceeds tens of EeV, and 

• neutrinos, with the first results from Superkamiokande, Baikal and AMANDA 
presented at this meeting. 
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Figure 2: Absorption of photons and protons on the interstellar light shown in Fig. 1. 
Shown is the absorption length in megaparsecs as a function of photon/proton energy. 

A novel, but essential aspect of the "New Astronomy" is that the Universe is not 
transparent to photons of TeV energy, and above. The same is true for protons once 
their energy exceeds 50 EeV. In the high energy sky, only neutrinos can reach us 
from the very edge of the Universe. The transparency of the Universe to photons 
and protons is shown in Fig. 2. Energetic photons are efficiently decelerated by pair 
production of electrons on background light above a threshold 



where E and e are the energy of the accelerated and background photon, respectively. 
Therefore TeV photons are absorbed on infrared light, PeV photons on the cosmic 
microwave background and EeV photons on radio-waves. It is, for instance, likely 
that absorption effects explain why Markarian 421, the closest blazar on the EGRET 
list at a distance of ~100 Mpc, produces the most prominent TeV signal. 
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Also protons interact with background light by photoproduction of the A- 
resonance just above the threshold for producing pions: 

2E p e > {ml - m 2 ) . (2) 

The major source of energy loss of ~100 EeV protons is photoproduction of pions 
on a target of cosmic microwave photons. The Universe is therefore opaque to the 
highest energy cosmic rays, with an absorption length of only tens of megaparsecs 
when their energy exceeds 10 20 eV. Lower energy protons, below threshold (0), do 
not suffer this fate. They can, however, not be used for astronomy because their 
directions are randomized in the microgauss magnetic field of our galaxy. 



2. Gravitational Waves |T[] 

The asymmetric collapse, e.g. of a rotating star, near the center of our galaxy will 
result into the supernova display astronomy is waiting for: the simultaneous obser- 
vation of light, neutrinos and gravitational waves could be the scientific event of all 
times. If we make the optimistic assumption that a similar amount of energy is 
emitted in gravitational waves and in light, i.e. one hundreds of a solar mass, gravi- 
tational antennas will detect a whopping signal of Sh = 10~ 18 . This deformation of 
the transverse components of the space-time tensor h^(x — ct) is detected at Earth 
in the form of gravitational waves. The passage of gravitational waves is revealed by 
a change in distance L between a pair of masses which are, ideally, separated by half 
the wavelength of the signal: 

5h= 6 4. (3) 

2 

Given the quadrupole structure of the waves, sensitivity is improved by using an 
orthogonal pair of spectrometers; see Fig. 3. 

Besides gravitational collapse, other anticipated sources of gravitational waves 
come in two categories: point sources and stochastic radiation. In the first category 
we identify the coalescence of a pair of neutron stars or black holes, the spin-down ra- 
diation of pulsars and binary systems. Examples of stochastic radiation include grav- 
itational waves emitted by early topological defects (today such radiation is strongly 
limited by the fact that it would destroy pulsar timing if present at significant lev- 
els) and the universal background of gravitons. Gravitons couple more weakly than 
neutrinos. In standard big bang cosmology they decouple at 10~ 22 seconds and, in 
contrast with the IK cosmic neutrinos, we can contemplate realistic methods for their 
detection. Any trace of gravitons would, of course, be wiped out during an epoch of 
inflation. 

Because interferometers operate by comparing candidate signals to templates of 
predicted signatures of gravitational waves, the theoretical computation of the defor- 
mation of the metric tensor is critical. Figure 4 displays the gravitational wave from 



4 



SCHEMATIC INTERFEROMETRIC DETECTOR 




PHOTODETECTOR 



LASER 



Figure 3: Schematic of LIGO. 




Figure 4: Gravitational wave from a binary. 



a binary system calculated to leading order in the potential 
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Such post-Newtonian calculations have been performed to order e 5 / 2 and results to 
third order will soon be available. 
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In the push to commission improved instruments we can identify three directions: 
i) coincident resonant bar detectors, ii) LIGO with a sensitivity in the 10~10 3 Hz 
frequency range, and iii) space-based interferometers such as VIRGO. Because time 
varying gravitational fields such as clouds, prevent the operation of ground-based 
detectors at frequencies below 1 Hz, space instruments are our only window to possible 
signals from supermassive objects such as the black holes which power active galaxies. 

LIGO construction has passed its two thirds mark towards completion, with 
progress measured in terms of dollars. With L = 4 km in Eq. (3), the instrument 
will initially reach a sensitivity of 5h > 10~ 21 . The challenge is to scale up the pro- 
totype, presently operating at Caltech, by a factor of 100 without loss of sensitivity. 
Optical precision has to reach one thousands of the micrometer wavelength of the 
laser light monitoring the distance between the pairs of masses in Fig. 3. Although 
a pair of neutron stars coalesces every minute, somewhere in the Universe, the signal 
of a nearby event is expected to be only of order 1 per year. The sensitivity of the 
instrument is limited by its frequency range to systems not much heavier than a solar 
mass. Yet black hole mergers may be viewed by LIGO out to 10~10 2 Mpc, at a rate 
of 10 per year, possibly more. 

3. Gamma Ray Astronomy on Earth and in Space |]2[] 

After two decades, ground-based gamma ray astronomy has become a mature science. 
The shower imaging method, pioneered by the Whipple telescope, emerged as the 
winning technique. Data taken on May 7, 1996 on the flaring active galaxy Markarian 
421 testify to this statement. Shown in Fig. 5 is the spectacular enhancement of 
photons in the direction of the source, over a small uniform background of sky light. 
At this meeting, appropriately, the Whipple group presented a string of interesting 
new results which include: 

• evidence that a source can be tracked to large zenith angle. This raises the 
threshold of the instrument and provides evidence that the Markarian 421 spec- 
trum extends well beyond 5 TeV, 

• evidence that this blazar emits TeV-photons in bursts with a duration of order 
a few days; see Fig. 6, 

• first hints of a correlation between the optical and TeV variability of blazar jets, 

• observation of a burst lasting only 15 minutes, suggesting emission from very 
localized regions of the galaxy, presumably the jet (more about jets later on), 

• identification of the first TeV source not detected in the GeV region by satellite- 
borne experiments such as EGRET: Markarian 501, another nearby blazar. 
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Figure 5: Signal and background of TeV photons during a flare of Markarian 421. 

Then there was the dog that didn't bark: no TeV emission from the supernova 
remnants IC433 and 7-Cygni, only upper limits from the Whipple and HEGRA detec- 
tors. If old supernova remnants are indeed the accelerators of the bulk of the cosmic 
rays, a healthy TeV signal is expected from the production and decay of neutral pions 
produced in the nuclear interactions of the nuclei with ambient matter in the shock. 
Such systems are a test-bed for models of shock acceleration and these observations 
did not go unnoticed. Models have already been tweaked in order to accommodate 
the upper limits in the TeV region. The conclusion is inevitably that accelerated 
electrons which scatter ambient light to GeV energy, are the origin of most of the 
flux observed by EGRET. Where do the cosmic rays come from? They are, most 
likely, only produced in the very late stages of the remnant, or in the rare supernova 
remnants with strong winds. 

The field of gamma ray astronomy is buzzing with activity to construct second- 
generation instruments. Space-based detectors are extending their reach from GeV 
to TeV energy with AMS and, especially, GLAST, while the ground-based Cherenkov 
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Figure 6: Time variation of the TeV photon flux from the blazar Markarian 421. 



telescopes are designing instruments with lower thresholds. In the not so far future 
both techniques should generate overlapping measurements in the 10 ~ 10 2 GeV 
energy range; see Fig. 7. While all ground-based experiments aim at lower threshold, 
better angular- and energy resolution, and a longer duty cycle, one can identify a 
three-prong attack on the construction of improved air Cherenkov telescopes: 

i. larger mirror area, exploiting the parasitic use of solar collectors during night- 
time (CELESTE and STACEY), 

it. better, or rather, ultimate imaging in the 17 m MAGIC mirror, 

in. larger field of view using multiple telescopes (VERITAS, HEGRA and TOKYO 
telescope arrays). 

VERITAS, for example, is an array of 9 upgraded Whipple telescopes, each with 
a field of view of 6 degrees. These can be operated in coincidence, or be pointed at 
9 different 6 degree bins in the night sky, thus achieving a very large field of view. 
The merits of stereo imaging, i.e. viewing the same air shower with multiple, imaging 
telescopes, was hotly debated. Encouraging initial measurements were presented by 
the TOKYO array, located in Utah, and by the HEGRA telescope array. The debate 
will be settled by experiment. 

There is a dark horse in this race: Milagro. The Milagro idea is to lower the 
threshold of conventional air shower arrays to 100 GeV by uniformly instrumenting 
an area of 10 3 m 2 or more (no sampling!). For time- varying signals, such as bursts, 
the threshold could be even lower. One can instrument a pond with photomultipliers 
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Figure 7: Thresholds of gamma ray telescopes. 

(Milagro), or cover a large area with resistive plate chambers (ARGO), or even with 
muon detectors (Hanul) which identify point sources of muons produced in photon- 
induced air showers. 

The conference was buzzing with rumors of observations of photons from active 
galaxies with energies as high as 50 TeV. Confirmation would revolutionize high 
energy astrophysics; we will return to this in the discussion of active galaxies. 

4. Proton Astronomy: EeV Cosmic Rays[§ 

Around 1930 Rossi and collaborators discovered that the bulk of the cosmic radia- 
tion is not made up of gamma rays. This marked the beginning of what was then 
called "the new astronomy", and we refer to as cosmic ray physics today. It is only 
"astronomy" above 5 x 10 19 eV or so, where the arrival directions of the charged 
cosmic rays are not scrambled by the ambient magnetic field of our own galaxy. As 
already mentioned, we suspect that the bulk of the cosmic rays are accelerated in the 
blastwaves of supernovae exploding into the interstellar medium. This mechanism 
has the potential to accelerate particles up to energies of 10 3 TeV where the cosmic 
ray spectrum suddenly steepens: the "knee" in the energy spectrum. We have no 
clue where and how cosmic rays with energies in excess of 10 3 TeV are accelerated. 
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We are not even sure whether they are protons or iron, or anything else. The origin 
of cosmic rays with energy beyond the "knee" is one of the oldest unresolved puzzles 
in science. 

It is sensible to assume that, in order to accelerate a proton to energy E, the 
size R of the accelerator must be larger than the gyroradius of the particle in the 
accelerating field B: 

E 

R > -Rgyro = • (5) 

This yields a maximum energy 

E < BR (6) 

by dimensional analysis and nothing more. The largest structures in our own galaxy 
are old supernova remnants with R ~ 10 2 parsecs. For the microgauss field of the 
galaxy B ~ 1CT 6 G, equation (5) yields a maximum energy of 10 5 TeV. Inefficiencies 
in the supernova shock acceleration mechanism will reduce this upper limit by a factor 
^shock/c — 0.1, as well as by another penalty factor for converting shock energy into 
acceleration, which is also of order 0.1. Therefore, the highest energy one can reach 
in practice is 10 3 TeV, i.e. the energy of the "knee" in the cosmic ray spectrum. 

Cosmic rays with energy in excess of 10 20 eV have been observed, some five orders 
of magnitude in energy above the supernova cutoff. To reach higher energy, one has to 
dramatically increase B and/or R. This argument is difficult to beat — it is basically 
dimensional. Although imaginative arguments actually do exist to avoid this impasse, 
it is generally believed that our galaxy is too small and its magnetic field too weak 
to accelerate the highest energy cosmic rays. Nearby active galactic nuclei distant 
by ~ 100 Mpc are obvious extra-galactic candidates. The jets of blazars support 
magnetic fields of 10 G over distances of 10~ 2 parsecs or more. Using Eq. (5) we 
reach energies of 10 20 eV, possibly higher because of beaming. Such speculations are 
reinforced by the fact that a cursory glance at the EGRET and Whipple results is 
sufficient to convince oneself that blazars are also the dominant (exclusive?) sources 
of the highest energy gamma rays. 

This raises the obvious question whether the highest energy cosmic rays point 
back to blazars? Astronomy with protons becomes possible once their energy has 
reached a value where their gyroradius in the microgauss galactic field exceeds the 
dimensions of the galaxy. From this point of view protons with 10 20 eV energy point 
at their sources with degree-accuracy. At this energy, their mean- free-path in the 
cosmic microwave background is unfortunately reduced to only tens of megaparsecs; 
see Fig. 2. A clear window of opportunity emerges: Are the directions of the cosmic 
rays with energy in excess of ~ 5 x 10 19 eV correlated to the nearest AGN (red-shift 
z less than 0.02), which are known to be clustered in the so-called "super-galactic" 
plane? Although far from conclusive, there is some evidence that such a correlation 
may exist, but not all experiments agree. Lack of statistics at the highest energies is 
a major problem. 
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Another problem is that the pointing accuracy is not really understood. It depends 
on the distance d to the source and the gyroradius in the intergalactic magnetic field: 
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Speculations on the strength of the inter-galactic magnetic field range from 10 to 
10~ 12 Gauss. For a distance of 100 Mpc, the resolution may therefore be anywhere 
from sub-degree to nonexistent. It is reasonable to expect that magnetic fields are 
higher in regions where matter is clustered. Higher values may therefore be ap- 
propriate for the local cluster and the super-galactic plane, regions relevant to this 
problem. Optimistically, one may anticipate that future high statistics experiment 
such as HIRES in Utah and the Auger giant air shower array will provide indirect 
information on the magnitude of the magnetic fields between galaxies. 

Before turning to experiment, a word about topological defects. Figure 8 shows the 
Fly's Eye cosmic ray spectrum in the EeV energy range, with the highest energy event 
ostensibly above a visual extrapolation of the lower energy events. One event has been 
sufficient to claim evidence for a top-down spectrum of "the" highest energy cosmic 
ray. It is suggested that the particle is the decay product of a 10 24 GeV (the GUT 
unification scale) topological defect such as a monopole, string.... This enthusiastic 
leap of faith can be excused. Topological structures are deeply connected to gauge 
theories and cannot be studied in accelerator experiments. Non-accelerator particle 
physics provides unique opportunities here. A topological defect will suffer a chain 
decay into GUT particles X,Y, which subsequently decay to the familiar weak bosons, 
leptons and quark-gluon jets. Cosmic ray protons are the fragmentation products of 
jets. We know from accelerator data that, among the fragmentation products of jets, 
neutral pions (decaying into photons) dominate protons by two orders of magnitude. 
Therefore, if topological defects are the origin of the highest energy cosmic rays, they 
must be photons. This is a problem because the highest energy event in Fig. 8 is not 
likely to be a photon. A photon of 3 x 10 20 eV interacts with the magnetic field of 
the earth far above the atmosphere and disintegrates into lower energy cascades; on 
average eight at this particular energy. The measured shower profile of the event does 
not support this assumption. One can live and die by a single event! 

Nevertheless consideration of topological defects deserves our attention because 
they can only be probed in non-accelerator experiments. That GUT theories are 
super symmetric (and so are half of the particles!) has been routinely forgotten in 
the literature. SUSY in the sky? Although supersymmetry modifies lots of factors 
in calculations that can, at best, be qualitatively tested, I actually do not know of 
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Figure 8: The highest energy cosmic ray flux measured by the Utah Fly's Eye. 

any clear SUSY signature in a cosmic ray detector. The stable, lightest LSP (the 
neutralino?) would be very much like a cosmic neutrino and leave no signature. 

Fig. 2 should remind us that, if the sources of cosmic rays are beyond 10 2 Mpc, 
conventional astronomy cannot find them. In this case, the absorption of the beam 
on the microwave background becomes a signature for distant sources, irrespective 
of the pointing precision of the cosmic rays at earth. The main problem today is 
statistics. After particles with energies in the vicinity of 100 EeV were discovered 
at Haverah Park, we accumulated a total of four events whose energy clearly exceed 
10 20 eV, using three different detectors: AGASA, Yakutsk and the Fly's Eye. The 
latter is being replaced by a technically superior instrument with larger collection 
area: the HIRES detector. Construction of a 10 4 km 2 array, one hundred times larger 
than the AGASA array operating in Japan, has been proposed and will hopefully be 
launched soon as the "Auger" project. 
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5. Intermezzo: Gamma Ray BurstsjQ 

The observation of gamma ray bursts confronts us with a seemingly straightforward 
problem: 

• a large amount of energy is released (of order 10 51 ergs, like a supernova), 

• in high energy photons (and only high energy photons — KeV and above), 

• in a very short time (typically seconds, or less). 

One burst a day is, on average, detected from a direction not correlated to our 
galaxy All evidence points at a cosmological origin of the sources, with an average 
redshift (z) ~ 1. 

Both the energetics and frequency point at neutron stars: a pair of neutron stars 
will merge once every minute, somewhere in the Universe. The rest is special relativity. 
High energy and short time can be accommodated if emission takes place in a site 
boosted towards the observer by a Lorentz factor 7. If the source moves towards us 
we expect that 

E = E ol n (9) 

and 

M = Rch n . (10) 

Here E Q is the energy in the frame of the source and R is its size. The power n 
is model-dependent and different for relativistic shocks in fireballs, jets, oscillating 
cosmic strings... For instance, models with n = 2 require 7 = 10 2 ~ 10 3 in order to 
describe the qualitative features of gamma ray bursts. 

As is the case for most acceleration problems in astrophysics, the devil is in the 
details. It has not been easy to produce the blueprint for a gamma ray burst. Al- 
though the details can be complex, the overall idea of fireball models is that a solar 
mass of energy is released in a compact region of radius R ~ 10 2 km. Only neutrinos 
escape because the fireball is opaque to photons. A significant fraction of the photons 
in the fireball is indeed above pair production threshold and produces electrons. It 
is straightforward to show that the optical depth of the fireball is of order 10 13 . It is 
theorized that a relativistic shock, with 7 ~ 10 2 or more, expands into the interstel- 
lar medium and photons escape only when the optical depth of the shock has been 
sufficiently reduced. There is, at least, one major problem: the shock is expected to 
lose its energy to the ubiquitous protons in the interstellar medium and fizzle long 
before it has the opportunity to produce the high energy gamma ray display. 

Personally, I find all models unconvincing because they do not naturally explain 
why the temporal structure of the spectrum of a gamma ray bursts is chaotic and to- 
tally different from burst to burst. At this meeting Dar presented a model overcoming 
this objection. 
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Dar postulates the formation of a pair of jets along the rotation axis of the merging 
neutron star pair, a configuration reminiscent of blazar jets. The jet material is 
presumably composed of protons and heavier nuclei, up to iron, just like cosmic rays. 
A high Z atom moving with a 7-f actor of order 10 3 will be photoexcited and ionized 
by the ambient starlight it encounters in the vicinity of the neutron star merger. From 
the atom's point of view, eV starlight is indeed boosted to KeV energy. KeV photons, 
emitted when the atoms relax, are detected as MeV gammas in the observer frame 
as a result of beaming. Thus two successive boosts turn starlight into MeV gammas 
according to Eq. (9), with n — 2. The baryons, which are a problem in conventional 
models, are now the solution! More importantly, because light is emitted whenever the 
shock encounters stars, the detected time profile of the burst represents a tomographic 
image of the galaxy in the vicinity of the neutron star merger. This explains the lack 
of organized structure; every burst is different. Monte Carlo simulations based on 
the distribution of stars near the center of our own galaxy produce bursts with a rich 
structure, indistinguishable from real data. Si non e vero, e ben trovato... Dar has 
not offered a detailed model for the propagation of the jet, the astrophysicists were 
skeptical. 

6. Intermezzo: Blazars|5|] 

EGRET and the air Cherenkov telescopes have put blazars at the focus of high en- 
ergy astronomy. EGRET has detected high energy gamma ray emission, in the range 
20 MeV-30 GeV, from over 100 sources. Of these sources 16 have been tentatively, 
and 42 solidly identified with radio counterparts. All belong to the "blazar" subclass, 
mostly Flat Spectrum Radio Quasars, while the rest are BL-Lac objects. In a unified 
scheme of AGN, they correspond to Radio Loud AGN viewed from a position illumi- 
nated by the cone of a relativistic jet. Moreover, of the five TeV gamma-ray emitters 
identified by the air Cherenkov technique, three are extra-galactic and are also nearby 
BL-Lac objects. The data therefore strongly suggests that the highest energy photons 
originate in jets beamed at the observer. Several of the sources observed by EGRET 
have shown strong variability, by a factor of 2 or so over a time scale of several days. 
Time variability is more spectacular at higher energies; see Fig. 6. 

Does pion photoproduction by accelerated protons play a central role in blazar 
jets? This question was hotly debated. If protons are accelerated along with electrons, 
they will acquire higher energies, reaching PeV-EeV energy because of reduced energy 
losses. High energy photons result from photoproduction of neutral pions by protons 
on the abundant UV photons in the jet. Accelerated protons thus initiate a cascade 
which dictates the features of the spectrum at lower energy. From a theorist's point 
of view the proton blazar has attractive features because protons, unlike electrons, 
efficiently transfer energy from the black hole in the presence of the high magnetic 
fields required to explain the confinement of the jets. The issue of proton acceleration 
can be settled experimentally because the proton blazar is a source of high energy 
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protons and neutrinos, not just gamma rays. 

First order Fermi acceleration offers a very attractive model for acceleration in 
jets, providing, on average, the right power and spectral shape. Confronted with the 
challenge of explaining a relatively flat multi-wavelength photon emission spectrum 
which extends to TeV energy, models have converged on the blazar blueprint shown in 
Fig. 9. Particles are accelerated by Fermi shocks in bunches of matter travelling along 
the jet with a bulk Lorentz factor of order 7 ~ 10. Ultra-relativistic beaming with 
this Lorentz factor provides the natural interpretation of the observed superluminal 
speeds of radio structures in the jet. In order to accommodate bursts lasting a day in 
the observer's frame, the bunch size must be of order TcAt ~ 10~ 2 parsecs. Here T is 
the doppler factor, which for observation angles close to the jet direction is of the same 
order as 7. These bunches are, in fact, more like sheets, thinner than the jet's width 
of roughly 1 parsec. The observed radiation at all wavelengths is produced by the 
interaction of the accelerated particles in the sheets with the ambient radiation in the 
AGN, which has a significant component concentrated in the so-called "UV-bump'.' 




black accretion disk 
hole 



Figure 9: Possible blueprint for the production of high energy photons and neutrinos 
near the super-massive black hole powering an AGN. Particles, accelerated in sheet 
like bunches moving along the jet, interact with photons, radiated by the accretion 
disk or, produced by the interaction of the accelerated particles with the magnetic 
field of the jet. 

In electron models the multi-wavelength spectrum consists of three components: 
synchrotron radiation produced by the electron beam on the I?-field in the jet, syn- 
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chrotron photons Compton scattered to high energy by the electron beam and, finally, 
UV photons Compton scattered by the electron beam to produce the highest energy 
photons in the spectrum. The picture has a variety of problems. In order to reproduce 
the observed high energy luminosity, the accelerating bunches have to be positioned 
very close to the black hole. The photon target density is otherwise insufficient for 
inverse Compton scattering to produce the observed flux. This is a balancing act, 
because the same dense target will efficiently absorb the high energy photons by 77 
collisions. The balance is difficult to arrange, especially in light of observations show- 
ing that the high energy photon flux extends beyond TeV energy. The natural cutoff 
occurs in the 10-100 GeV region. Finally, in order to prevent the electrons from 
losing too much energy before producing the high energy photons, the magnetic field 
in the jet has to be artificially adjusted to less than 10% of what is expected from 
equipartition with the radiation density. 

For these, and other reasons already mentioned, the proton blazar has been devel- 
oped. In this model protons as well as electrons are accelerated. Because of reduced 
energy loss, protons can produce the high energy radiation further from the black hole. 
The more favorable production-absorption balance far from the black hole makes it 
relatively easy to extend the high energy photon spectrum above 10 TeV energy, even 
with bulk Lorentz factors that are significantly smaller than in the inverse Compton 
models. Because the seed density of photons is still much higher than that of target 
protons, the high energy cascade is initiated by the photoproduction of neutral pions 
by accelerated protons on ambient light via the A resonance. 

Model-independent evidence that AGN are indeed cosmic proton accelerators can 
be obtained by observing high energy neutrinos from the decay of charged pions, pho- 
toproduced along with the neutral ones. The expected neutrino flux can be estimated 
in six easy steps. 

i. The size of the accelerator R is determined by the duration, of order 1 day, over 
which the high energy radiation is emitted: 

R = Ftc = 10~ 2 parsecs for t = 1 day and T = 10. (11) 

ii. The magnitude of the S-field can be calculated from equipartition with the 
electrons, whose energy density is measured experimentally: 

B 2 

- — = p(electrons) ~ 1 erg/cm 3 . (12) 
2/i 

This yields a value for the magnetic field of 5 Gauss. 

Hi. In shock acceleration the gain in energy occurs gradually as a particle near the 
shock scatters back and forth across the front gaining energy with each transit. 
The proton energy is limited by the lifetime of the accelerator and the maximum 
size of the emitting region, R 

E < KZeBRc . (13) 
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Here Ze is the charge of the particle being accelerated and B the ambient 
magnetic field. The upper limit basically follows from dimensional analysis. 
It can also be derived from the simple requirement that the gyroradius of the 
accelerated particles must be contained within the accelerating region R; see 
Eq. (5). The numerical constant K ~ 0.1 depends on the details of diffusion 
in the vicinity of the shock, which determine the efficiency by which power in 
the shock is converted into acceleration of particles. In some cases it can reach 
values close to 1. The maximum energy reached is 

£max = eBRc = 5 x 10 19 eV 

for B = 5 Gauss and R = 0.02 parsecs. We here assumed that the boost of 
the energy in the observer's frame approximately compensates for the efficiency 
factor, i.e. KT ~ 1. 

The neutrino energy is lower by two factors which take into account i) the 
average momentum carried by the secondary pions relative to the parent proton 
((xp) ~ 0.2) and ii) the average energy carried by the neutrino in the decay 
chain n + — > — > e + u e u fl1 which is roughly 1/4 of the pion energy because 
equal amounts of energy are carried by each lepton. The maximum neutrino 
energy is 

£,max = E max (x F ) i ~ 10 18 eV , (14) 

i.e. neutrinos reach energies of 10 3 PeV. 

iv. The neutrino spectrum can now be calculated from the observed gamma ray 
luminosity. We recall that approximately equal amounts of energy are carried 
by the four leptons that result from the decay chain n + — > — > e + v e V lx . 
In addition the cross sections for the processes jtry — > pit and £>7 nir + at 
the A resonance are in the approximate ratio of 2:1. Thus 3/4 of the energy 
lost to photoproduction ends up in the electromagnetic cascade and 1/4 goes 
to neutrinos, which corresponds to a ratio of neutrino to gamma luminosities 
(L„/L 7 ) of 1/3. This ratio is somewhat reduced when taking into account that 
some of the energy of the accelerated protons is lost to direct pair production 
(p + 7 — > e + e~p): 

L„ = Al 7 . (15) 

In order to convert above relation into a neutrino spectrum we have to fix the 
spectral index. We will assume that the target photon density spectrum is 
described by a E~( 1+a ^> power law, where a is small for AGN with flat spectra. 
The number of target photons above photoproduction threshold grows when 
the proton energy E p is increased. If the protons are accelerated to a power 
law spectrum with spectral index 7 = 2, the threshold effect implies that the 
spectral index of the secondary neutrino flux is also a power law, but with 
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an index flattened by (1 + a) as a result of the increase in target photons at 
resonance when the proton energy is increased: 



dN v 
dE„ 



E 

max 



(1-a) 



(16) 



For a flat photon target with a — 0, the neutrino spectrum will flatten by just 
one unit giving E^f- ~ constant. From Eqs. (15) and (16) 

£„ ma x dN„ E 2 3 



/ 



u Assuming that the high energy 7 ray flux from Markarian 421 results from 
cascading of the gamma ray luminosity produced by Fermi accelerated pro- 
tons, we obtain the neutrino flux from the measured value of L 1 of 2 x 
1(T 10 TeV cm' 2 s -1 : 



dN u 3 L 7 1 + a 



dE v 1 3 E v max E u 



E v 



E 

max 



5 x lO-^cm^s- 1 , . 
E ' (18) 



where the numerical estimate corresponds to a = and the value of E v max of 
Eq. (14). The neutrino flux is essentially determined by the value for E v max . 

vi. In order to calculate the diffuse flux from the observed blazar distribution, 
we note that the EGRET collaboration has constructed a luminosity function 
covering the observation of the ~20 most energetic blazars and estimated the 
diffuse gamma ray luminosity. From the ratio of the diffuse gamma ray flux 
and the flux of Markarian 421, we obtain that the effective number of blazars 
with Markarian 421 flux is ~130 sr _1 . The diffuse neutrino flux is now simply 
estimated by multiplying the calculated flux for Markarian 421 by this factor. 

This concludes our calculation. It illustrates how the proton, unlike the electron, 
blazar requires no artificially large T factors and no fine-tuning of parameters. For 
the proton blazar, radiation and magnetic fields are in equipartition, the maximum 
energy matches the BR value expected from dimensional analysis and, finally, the 
size of the bunches is similar to the gyroradius of the highest energy protons. It is 
not a challenge to increase gamma ray energies well beyond the TeV energy range. 
Reasonable variations of the values of magnetic field strength B, the efficiency pa- 
rameter K and the Doppler boost factor T may allow us to account for the highest 
energy cosmic rays with E ~ 3 x 10 20 eV. 

The probability to detect a TeV neutrino is roughly 1CT 6 . It is easily computed 
from the requirement that, in order to be detected, the neutrino has to interact within 
a distance of the detector which is shorter than the range of the muon it has produced. 
Therefore, 

P^^^czAEZ, (19) 

^int 
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where i? M is the muon range and A; nt the neutrino interaction length. For energies 
below 1 TeV, where both the range and cross section depend linearly on energy, n = 2. 
Between TeV and PeV energies n = 0.8 and A = 10~ 6 , with E in TeV units. For EeV 
energies n = 0.47, A = 10~ 2 with E in EeV. 

We are now ready to compute the diffuse neutrino event rate by folding the neu- 
trino spectrum of Eq. (18) with the detection probability of Eq. (19). We also multiply 
by 130 sr -1 for the effective number of sources: 

/•El/max (] /V „ 

<f>" = J -^rP v ^(E v )dE v ~ 40 km" 2 year" 1 sr" 1 . (20) 

which implies a yield of two neutrinos every three days in a kilometer-scale detector, 
assuming only 27r coverage. 

We have already drawn attention to the 10 TeV maximum photon energy as the 
demarkation line between the electron and proton blazars. The ~10 GeV cutoff in 
the inverse Compton model can be pushed to the TeV range in order to accommodate 
the Whipple data on Markarian 421, but not beyond. Bringing the accelerator closer 
to the black hole may yield photons in excess of 10 TeV energy — they have, however, 
no chance of escaping without energy loss on the dense infrared background at the 
acceleration site. 

HEGRA has been monitoring the 10 closest blazars, including Markarian 421, 
with its dual telescope systems: the scintillator and naked photomultiplier detector 
arrays. They announced at a variety of meetings (but not at this one!) that their 
upper limit on the flux of photons of 50 TeV energy and above, for the aggregate 
emission from the ten nearest blazars, may actually be a signal. This could provide 
the first compelling evidence that blazar jets are indeed true proton accelerators. 



7. The Birth of High Energy Neutrino Astronomy [|| 

With the rapidly expanding Baikal and AMANDA detectors producing their first hints 
of neutrino events, observation of neutrinos from AGN could establish the production 
of pions and identify the accelerated proton beams as the origin of the highest energy 
photons and the highest energy cosmic rays. A definite answer may not be known 
before neutrino telescopes reach kilometer size. The builders should take note that, 
although smaller neutrino fluxes are predicted than in the generic AGN models of a 
few years ago, they are all near PeV energy where the detection efficiency is increased 
and the atmospheric neutrino background negligible. Because of the beaming of the 
jets, the neutrinos have a flat spectrum peaking near the 10 6 TeV maximum energy. 
The actual event rates are, in the end, not very different. 

These models strongly favor the construction of neutrino telescopes following a 
distributed architecture, with large spacings of the optical modules and relatively 
high threshold: ~ 100 meter spacing among strings and ~ 20 meter spacings of 
OMs along a string (OM: pressure vessel containing a conventional photomultiplier 
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tube and, possibly, data acquisition electronics). Models where most neutrinos have 
very high energy, open up opportunities for alternative techniques such as the radio 
technique, or the detection of horizontal air showers with giant air shower arrays. Both 
were extensively discussed at this meeting. Optimists, on the other hand, can find 
reasons to anticipate the discovery of AGN neutrinos with much smaller telescopes, 
such as the existing Macro, Baikal and AMANDA detectors. With a sufficiently high 
proton target density in the acceleration region, much larger fluxes of neutrinos may 
be produced in a proton-proton cascade. The predicted fluxes are however model- 
dependent. It is also possible, even likely, that photons do not escape the source or, 
escape after significant energyloss. Such absorption effects increase the neutrino flux 
relative to the observed high energy photon flux, also leading to the prediction of 
larger neutrino fluxes. 

The motivation for building neutrino telescopes is not limited to the observation 
of AGN. More than 6 orders of magnitude in photon energy, or wavelength, are left 
unexplored between GeV photons and EeV protons in Fig. 1. Neutrino telescopes 
have been conceived to fill this gap. Given the history of astronomy, it is difficult 
to imagine that this will be done without making major, and most likely totally 
surprising, discoveries. The 18 orders of magnitude in wavelength, from radio-waves 
to GeV gamma rays, are indeed sprinkled with unexpected discoveries. Neutrinos 
have the further advantage that, unlike photons of TeV energy and beyond, they are 
not absorbed by interstellar light. They can, in principle, reach us from the edge of the 
Universe. If, for instance, the sources of the high energy cosmic rays are well beyond 
the Virgo cluster, the photon window for their exploration closes above 100 TeV. 

Detectors may have to reach kilometer scale before finding the accelerators of the 
highest energy cosmic rays, improving significantly the searches for cold dark matter 
particles or WIMPs, and to search meaningfully for the cosmic sources of gamma 
ray bursts. In order to explore the plausible region of astronomical parameter space 
for these fascinating cosmic enigmas, a high energy neutrino telescope must contain 
several thousand optical modules in a volume of order 1 kilometer on a side. Model 
building suggests that some, and most likely the most exciting, discoveries may be 
within reach of much smaller detectors with effective telescope area of order 10 4 m 2 . 
Experience with small detectors is, in any case, an important intermediate step and 
early indications are that a kilometer-scale detector could be constructed in 5 years 
using existing technologies. 

Consisting of several thousand OMs deployed in natural water or ice, even the 
ultimate scope of these detectors is similar to that of the Superkamiokande exper- 
iment, which presented confirmation of the Kami oka solar neutrino results at this 
meeting. Being optimized for large effective area rather than low threshold (GeV 
or more, rather than MeV), they are complementary to SuperK. The challenge to 
deploy the components in an unfriendly environment is, however, considerable. With 
a price tag which may be as low as a relatively cheap fixed-target experiment at an 
accelerator, but could be as high as that of a LHC detector, this must be one of the 
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best motivated large-scale scientific endeavors ever. 

As for conventional telescopes, at least two are required to cover the sky. As 
with particle physics collider experiments, it is very advantageous to explore a new 
frontier with two or more instruments, preferably using different techniques. This 
goal may be achieved by exploiting the parallel efforts to use natural water and ice 
as the Cherenkov medium for particle detection. 

There has been a heightened level of activity in this field in the last 3 years. In 
the Spring of 1993, using the frozen ICG clS du platform for easy deployment, the Lake 
Baikal group deployed a small telescope consisting of 36 optical modules. They plan to 
complete the detector consisting of 200 optical modules by 1998. The Russian-German 
collaboration presented the neutrino-induced, upgoing muon event shown in Fig. 10. 
With 19 channels reporting on 4 strings, it is gold-plated. Deep ocean water should 
be superior in optical quality to that in Lake Baikal. Two collaborations, ANTARES 
(France) and NESTOR (Greece), are developing the infrastructure and technologies 
for the deployment of neutrino telescopes in the Mediterranean basin. It is, at this 
point, only fitting to recall the pioneering role played by the now defunct DUMAND 
experiment in Hawaii. They made key conceptual and technological contributions to 
this field. 

Using natural Antarctic ice as a particle detector, the AMANDA A detector of 
80 OMs, positioned near 1 kilometer depth, has been taking data for almost 3 years. 
A deep array of ~ 300 OMs, AMANDA B, was completed a few days after this 
conference. In the next season the first strings of kilometer length will be deployed, 
with the goal to commission a kilometer-scale detector over a period of approximately 
5 years. The problem of reconstructing muons in a kilometer-scale detector has been 
assessed experimentally by i) studying muon tracks registered in both the 1 and 2 
kilometer-deep detectors and ii) investigating linear energy response in AMANDA A. 

Coincidences between AMANDA A and B are triggered at a rate of 0.1 Hz. Every 
10 seconds a muon is tracked over 1.2 kilometer; a typical event is shown in Fig. 11. 
Below 1500 m the vertical muon triggers 2 strings separated by 79.5 m. The distance 
along the Cherenkov cone is over 100 m, yet, despite some evidence of scattering, 
propagation of the muon at the speed of light can be readily identified. This, and 
several other calibration methods indicate a scattering length 2 orders of magnitude 
larger in the deep detector than in the shallow one. The collaboration has analysed 
over 5 x 10 5 A-B coincidences with no evidence for any source of misreconstructed 
background events. Reconstruction of such tracks to a degree in the bubble-free ice 
should not represent a challenge. 

Remnant bubbles, much larger in size than expected, limit the scattering length 
to tens of centimeters in the shallow detector, preventing the reconstruction of muon 
tracks. With this set-back nature did however provide the AMANDA collaboration 
with the hint that scattering can be exploited to measure energy. Bubbles very 
effectively diffuse and contain the light from showers, e.g. initiated by electron neu- 
trinos. The shower can be mapped and its energy measured provided the point of 
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Figure 10: An upgoing muon initiated by a neutrino interaction below the Lake Baikal 
detector. The Cherenkov signal is mapped by 19 optical modules. 

origin is within 50 m of the instrumented ice. Candidate atmospheric neutrino events 
have been identified in the range of 100 GeV to 1 PeV. A measurement of the spec- 
trum requires a detailed understanding of all systematics. AMANDA A is the ideal 
prototype to develop the methods for measuring energy. It is, with its reduced de- 
tector dimension and scattering length, a scale model of a kilometer-cubed detector 
in bubble-free ice. 

Let me conclude by trying to infuse some sanity in the non-debate on "water and 
ice". It is a non-debate because, ideally, we want both. Given the pioneering and 
exploratory nature of the research, we most likely need both. All indications are 
that water and ice have complementary optical properties: while the "attenuation" 
lengths are comparable for the blue wavelength photons relevant to the experiments, 
attenuation is dominated by scattering in ice and by absorption in water. Both have 
a problem: scattering in ice, potassium decay and bioluminescence in water. Both 
problems can be solved. The high rates of noise events, especially in a kilometer-size 
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Figure 11a: Cosmic ray muon track triggered by both AMANDA A and B. Trigger 
times of the optical modules are shown as a function of depth. The diagram shows 
the diffusion of the track by bubbles above 1 km depth. Early and late hits, not 
associated with the track, are photomultiplier noise. 
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Figure lib: Cosmic ray muon track triggered by both AMANDA A and B. Trigger 
times are shown separately for each string in the deep detector. In this event the 
muon mostly triggers OMs on strings 1 and 4 which are separated by 79.5 m. 
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water detector, can be removed by only triggering on pairs of OMs. The solution to 
the scattering problem in ice has been demonstrated by Monte Carlo simulation. It 
will be tested on the data which are collected at a rate of ~100 Hz while I am writing 
this. Because of the fantastic transparency of ice, ~20 OMs report in a typical event, 
compared to 6~10 in a water detector with the same effective area. The first photon 
reaching 5, sometimes more, of the 20 OMs will not be scattered. Their arrival times 
and the information from the other OMs are merged in a likelihood fit which delivers 
track reconstruction to better than 2 degrees. 

A proposal for a surface detector Hanul, which means "sky" in Korean, introduces 
new technology in this neutrino race. It will be composed of modules, somewhat larger 
than Macro and LVD, which measure muon charge and energy with a large magnet 
sandwiched between rpc's and Cherenkov counters. 

With SuperK, Macro, AMANDA and Baikal this promises to be a happy z/-year. 
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